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An iterative technique is developed for retrieval of 
thicknesses of selected atmospheric layers from VTPR "clear- 
column" radiance measurements. Layer mean temperatures for 
a simplified atmospheric model are retrieved by direct solu- 
tion of the radiative transfer equation, and are then used 
to compute thicknesses of key atmospheric layers bounded by 
commonly used pressure levels. The retrieval technique 
‘illustrates the use of reference wave numbers that vary from 
layer to layer, Transmittance tuning is employed to correct 
syatenatic errors in the retrieved mean temperatures. 
Thicknesses of key layers retrieved by the technique from 
"clear-column" radiances observed during a 24 hour period at 
scan spots between 15 N and 45 N are separated into three _ 
latitude-band samples. Each sample is subjected to stepwise 
multiple regression analysis to determine the thickness- 
specification of various standard layers in terms of the 
clear column radiances. RMS error-analyses resulting from 
the regression are then used to determine the quality of 
thickness-specifications of simulated tropospheres and 


stratospheres. 
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I. INTRODUCTION 


Retrieval of meteorological parameters from satellite 
radiance measurements has been the object of numerous 
research studies since Kaplan [1959] demonstrated that 
vertical temperature profiles of the atmosphere could be 
inferred from satellite spectral radiance observations in 
the 15 um band of carbon dioxide emission. Application of 
retrieval techniques to satellite measurements has been 
possible since the launch of NIMBUS III in 1969 and to date 
emphasis has been on retrieval of temperature profiles. 

However, for purposes of numerical weather prediction, 
the atmospheric thickness of a specified pressure interval 
is normally a more basic parameter than temperature (Fleming, 
1972]. Furthermore, thicknesses of specific atmospheric 
pressure intervals, or layers, may be retrieved by direct 
solution of the radiative transfer equation, thereby 
eliminating the need for "a priori" statistical information 
necessary for retrieval by regression or inverse matrix 
methods [Fritz, Wark, et al., 1972]. 

A direct retrieval method for obtaining specific layer 
mean temperatures, and hence thicknesses, from satellite 
radiance measurements was presented by Fleming [1972]. The 
method proved to be too time-consuming to be operationally 
adaptable by the National Environmental Satellite Service 


(NESS) to the National Meteorological Center (NMC) analysis 
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scheme as described by McMillin, Wark, et al. [1973], and 


the retrieved thickness values showed little improvement in 
accuracy when compared to results obtained by a statistical 
regression technique which gave a T(p) profile at mandatory 
levels from which thicknesses were computed [Smith and 
Fleming, 1972]. However, the concept proposed by Fleming 
[1972] appeared to have at least diagnostic merit for 
application to vertical analysis schemes such as that of 
Fleet Numerical Weather Central (FNWC), which employs a 
layer thickness as an input parameter [Holl, et al., 1964]. 
Fleming proposed a method whereby the unit square wave 
function denoted by R(ln p/P.) could be determined as a 
linear combination of spectral radiances, I,; with coeffi- 


cients, Cy chosen to minimize the right side of 


N Ps 
Ee, = AT Sf R[1n p/p,; Py2PoJd(in p) . a) 
11 p 
o 
Here 
AT = mean temperature for pressure interval (Py »Po)- 


Po» Ps = pressures at top of atmosphere and surface of 
earth. 


i=i1, 2, ..., N are the spectral intervals or 
channels. 
The R-functions were sought to provide, in sequence, exact 
fits to the square wave in the significant layers of the 
atmosphere. By Fleming's theory they should ideally fit 
the square wave of Fig. 1 between Py and Po and be zero 


elsewhere along the p-axis. 


cm! 





FIG. 1. Exact vs. calculated rectangular function 
R. Desired solution indicated by dashed line, 
calculated solution indicated by solid line. 


In practice, the calculated R functions did not fit the 
square wave exactly and left errors in adjacent pressure 
regions as depicted by areas A and B of Fig. 1. Inexactness 
of fit at the top of the square wave resulting in area C was 
considered unimportant because it tended to average out to 
the required thickness constraint over (Py »Po)- 

The implication was that the better the R function fits 
the square wave for a given layer (Py oPo)s the better the 

mean temperature of the layer was specified by the radiance 


a measurements. 
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The object of this thesis was first to develop a direct 


retrieval Sachudaus for layer-mean temperatures (and there- 
fore thicknesses) which required a minimum of computer time, 
and could conceivably be less sensitive to channel noise 
than Fleming's method. A second objective was designed to 
test retrieved thicknesses for optimum layer specification 
by statistical methods, as implied by Fleming's square wave 
concept. 

In this thesis a new simplified iterative retrieval 
technique was developed for use with VTPR "clear-column" 
radiances to obtain mean temperatures over greater pressure 
intervals than the mandatory pressure increments used by 
Fleming. An iterative retrieval technique was used because 
of the minimum computer time and space required. The use 
of larger pressure intervals was employed to reduce the 
effect of channel noise and to improve accuracy, since the 
accuracy of mean temperatures derived from radiance measure- 
ments usually improves with increasing pressure interval 
(Hayden, 1971]. 

Key-layer thicknesses derived hypsometrically from the 
mean T(p) profile were then subjected to stepwise multiple 
regression analysis to determine which key Tagen was best 
specified by the "clear-column" radiances. The troposphere 
and lower stratosphere were considered separately, that is 
layers crossing the tropopause level such as 500 to 50 mb 


‘were not considered. 
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Results of the regression analysis were examined first 


to determine which particular layer in the troposphere (and 
also which layer in the lower stratosphere) could be used 
as the most effective thickness from VIPR retrieval for use 
in the FNWC vertical-structure analysis scheme, which 
currently makes use of the 1000 to 300 mb thickness in 
processing of conventional sounding data. 

Next, combinations of sub-layers which span the tropo- 
sphere were examined to determine if a better statistical 
fit for the retrieved tropospheric thickness could be 
obtained from the combinations of sub-layers, and if so, 
then which combination was best. A similar study was made 


for the lower stratosphere. 
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II. SATELLITE DATA 


The launch of NOAA-II with its Vertical Temperature 
Profile Radiometer (VTPR) instrument in October, 1972 
marked a major improvement in radiance measurement capability. 
The VTPR instrument is superior to its earlier counterparts, 
the SIRS-A of NIMBUS III and the SIRS-B of NIMBUS IV, in that 
it has much better spatial resolution. Sub-satellite dimen- 
sions for scans spots are approximately 69 by 67 km for VTPR 
compared to 225 by 225 km for SIRS-A and SIRS-B. The 
improved resolution permits more accurate computation of 
"clear-column" radiances, which are equivalent to radiances 
that would be observed in completely clear skies. Use of 
"clear-column" radiances for retrieval eliminates the neces- 
sity to correct for cloud cover which was a significant 
problem in retrieval techniques developed for use with SIRS-A 
and SIRS-B data. 

NOAA~II orbits the earth every 115 minutes at an altitude 
of 1464 km. Figure 2 illustrates the earth projection of 
seven orbits. North to south portions of the orbit are 
indicated by solid lines; south to north portions are 
indicated by dashed lines. Shaded areas depict areal 
coverage during two orbits. North to south equator crossings 
occur at 0900 and south to north crossings occur at 2100 local 
solar time. 

The VTPR instrument scans perpendicular to the satellite 


path in 23 discrete steps from left to right, representing 


15 
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Satellite tracks for NOAA II VIPR coverage 








a scan path of 30.3 degrees both sides of the local nadir. 


Each step provides a "scan spot" which is observed by the 

VIPR instrument in six spectral intervals of the 15 um band of 
carbon dioxide, inone interval of the 12 um atmospheric window, 
and in one interval of the 19 um water vapor band. 3 

The elght spectral intervals, or channels, along with their 
respective half-widths and central wave numbers are listed 


in Table 1. 


TABLE 1. VTIPR channel designators corresponding half-widths, 
and central wave numbers. , . 


Channel 1 2 3 4 5 6 7 8 
ea 668.5 677.5 695.0 708.0 725.0 747.0 535.0 835.0 
cm 


Half-width 3.5 10.0 10.0 10.0 10.0 10.0 18.0 10.0 
(cem-1) 


Channels one through six are carbon dioxide channels, 
channel seven is the water vapor channel, and channel eight 
is the window channel. fies add relative measurement error 
between any two channels except the 668.5 em 72 channel is 
0.25 mwW/(m* ster em 1); maximum relative error between the 
668.5 com + channel and any other channel is 0.75 mW/(m? ster 
em~-), 

Other instruments aboard NOAA-II include a two-channel 
‘Scanning Radiometer (SR) which measures radiances in the 


10.4-12.5 and 0.5-0.7 um intervals. Resolution is more 
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refined in the SR scan spots than that of the VTPR, with sub- 


satellite dimensions at the nadir being approximately 7.5 by 
7.5 km. Statistical techniques are used to identify scanning 
radiometer measurements which signify cloud-free areas, and 
these "clear—column" SR window channel radiances are then 
used to determine sea surface temperatures for the VTPR 

Scan spots. 

Conversion of the VIPR raw radiance measurements to 
"clear-column" radiances is accomplished by first dividing 
Scan spots into analysis arrays. Scan spots from eight 
Successive scan lines are divided into three boxes, or sub- 
arrays, of 8 by 8, 8 by 7, and 8 by 8 spots as illustrated 
in Fig. 3, which is an enlargement of the boxed area outlined 
in Fig. 2. From the VIPR raw radiance measurements and the 
SR derived sea surface temperatures of scan spots within each 
Ssub-array, a single set of "clear-column" radiances is com- 
puted by statistical methods [McMillin, Wark, et al, 1973], 
and assigned to central scan spot locations indicated by 
X's in Fig. 3. 

The procedure is essentially to compute an 835.0 emt 
window channel radiance by the Planck formula using the sea 
surface temperature, and to compare the computed value 
against the measured 835.0 em? window channel radiance. 

If the measured 835.0 em! radiance value equals or exceeds 
the computed value (from the known sea-surface temperature 
field, denoted by SST), the radiances are considered to be 


in agreement and the scan spot is assumed to be cloud-free. 
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FIG. 3. VTPR scan pattern and data analysis array for 
the outlined box in Fig. e. 


If the computed value (from the SST) exceeds the measured 
value, the scan spot is assumed to contain significant 
cloud cover and equivalent "clear-column" radiance values 


are calculated by noting that for a given wave number, V5» 


I (v ) -I,(,) I (v,) -T 04) 


clr‘ ‘i _ elr F (2) 
Torpedg) ~2atvg) — Toy,0¥g) 1, (vg) 
toi V8 14 V8 Toir V8 qy Vg 
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Toi (vg) = window radiance computed from the SR 
r derived sea surface temperature. 


I, (v4 )315(v,) = raw radiances in channel 1 measured at 
scan spots 1 and 2 having the same sea 
surface temperature. 

Tein’) = desired "clear-column" radiance value 
for wave number 1° 

I, (vg) 15 (vg) = radiance measured for the 835.0 em 
window channel at scan spots 1] and 2. 


As shown in Fig. 4, the three points (I yp (¥,) 27 (vg)]; 


(I, (v,),T,(vg)], and [I,(v,),I5(vg)] lie on the straight line 


clr 


determined from the measured values of I,(v5)5 I, (vg), In(v4)5 
and I5(vg). The slope of this line is expressed by the right 
side of eq. (2). 

The value of Top’) can then be determined from the 
known value of Toirl’g): When computing Tip Ys? values, 
radiances from adjacent scan spots with different nadir 
angles are adjusted to a common zenith angle. 

After clear radiance values are obtained for all channels 
at each scan spot within a sub-array, all values from adjacent 
scan spots are examined. The sub-array maximum radiance set 
yields a single set of eight "clear-column" radiances, and a 
sea surface temperature positioned at the center of the sub- 
array. This results in reduction of raw radiances for the 
184 scan spots of the 8 by 23 analysis array, to "clear- 
column" radiances at only the centers of the three sub-arrays. 

The resulting VITPR "clear-column" radiance values and 
“associated data are recorded on archival tapes. An archival 


tape for April 12, 1973, was provided for this study through 


20 






(OVERCAST) 





es ke is aes oa ee ae es aes lees ee I es 








J (Y,) —— lor (Vg) 





FIG. 4. Procedure for determining "clear~column" 
radiances. 


the kind auspices of Dr. D. Q. Wark of NESS. From a printout 
of archival file II (the "clear column" file), radiances for 
the six carbon dioxide channels and the sea surface tempera- 
ture, as well as the geographic coordinates, were extracted 
for clear column scan spots between 15 and 45 degrees north 
latitude. 

Data from the tape printout were then transformed to 


correct dimensions as follows: 


N. Latitude in degrees = (tape value) x .10 -90 
W. Longitude in degrees = (tape value) x .10 
Radiances in mW/ (m® ster om7+) = (tape value) x .05 

, Sea surface temperature in °K = (tape value) x .20 + 269.9 
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Additional information on characteristics of the NOAA-II 


Satellite and its VTPR and SR radiance data may be found in 
NOAA technical reports [Fritz, Wark, et al., 1972] and 
[McMillin, Wark, et al., 1973] from which most of the details 


of this section were taken. 
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III. RETRIEVAL TECHNIQUE 


Given VTPR "clear-column" radiance values for channels 
one through six and the sea surface temperature computed 
essentially from SR measurements, corresponding layer mean 
temperature profiles were obtained using an iterative tech- 
nique derived from that initially presented by Smith [1970] 
and modified by Martin [1973] for the purposes of this 
study. Mean temperatures are then converted to layer 
thicknesses by use of the hypsometric equation, as if the 
mean temperatures of the iterative procedure applied to 


their central pressure levels. 


A. MATHEMATICAL DEVELOPMENT 

For a cloudless, non-scattering atmosphere in thermo- 
dynamic equilibrium, the spectral radiance observed at the 
top of the atmosphere for each channel is related to the 
vertical temperature profile and absorbing gas structure by 


the radiative transfer equation (RTE) [Fritz, Wark, et al., 


1972): 
x(po) dt, (p) 
I, = B, [T(p,) It, (pg) 5 3 oe B,[T(p) J ax(py ox?) (3) 
s 


(A) (B) 
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I, = spectral radiance in channel 1, (LET 23600560) 


B, (T(p) J 


Planck radiance function for channel i and 
temperature T at pressure level p, 


tT, (p) = fractional transmittance of the atmospheric 
CO, in channel i from pressure level p to 


Po = .O1 mb, 


x(p) 
in the vertical similar to log P/Po: 


Term (A) is the atmospheric transmittance of the Planckian 


radiance from the surface of the earth. Term (B) is the 


an arbitrary function of pressure which behaves 


atmospheric contribution to the radiance. Subscripts s and 


o refer to surface of the earth (1000 mb) and top of the 
atmosphere (.01 mb), respectively. 
When 100 pressure levels are linearly scaled by pe/T 


and the result adopted for x(p), it follows that 
p(J) = .01[1 + (J-1)(0.26087836)]!/2 (4) 


and J = 1,2,...,100 are pressure levels numbered from top 
of the atmosphere to surface of the earth, the radiative 
transfer equation can be rewritten [Martin, 1973] 


J=01 dt, (J) 


I, = B,[T(100)]1, (100) + : B,[T(5)] —33—— ad (5) 


f 
i =100 


The Planck radiance function is defined as 
3 (Cov4/T(S)) , 
B,(T(3)] = yyy te - 1) (6) 
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Vv, = wave number for channel i, 


1.9061 x 107° erg cm* sect ster, 


Qa a 
t ul 


= 1.43868 cm °K. 


When the 100 pressure levels are combined into the 17 
atmospheric layers depicted in Fig. 5, eq. (5) can be 


evaluated in quadrature form 


17 
i, = B, [T(100) Jt, (100) + a B, (K) At, (K) - (7) 
where 
By (i) = ~ {B,[T(01)] + 4B, [T(02.5)] + ByLT(O4) I} (8) 


is the layer mean Planck function for layer K = 1, and 
— a © 
By(K) = ¢ {B,[T(Jy-3)] + 4B (T(J,) + BLT(Sy+3)]} (9) 
is the layer-mean Planck function for layer K > 1. In (7), 
At, (K) is the transmittance in channel i from layer K as 
defined by 

At, (K) = Ts (Jy-3) = Ts (Jy+3) e (10) 


Finally, J, is the value of J at the center of layer kK. 


k 
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1000 MB 


100 


FIG. 5. Seventeen K-layer atmosphere model used 
for mean-temperature retrieval. Note that the 
top layer spans only four J-levels while all : 
others span seven J-levels. 





a Equation (7) can be written in iterative form as 


- I, - Tf = {By [1(100)] - BpLt(100) J}r, (100) 
(11) 
ar BML (x) - By(K)}At, (K) 
K=1 
where 
I, = observed radiance for channel i, considered to 


be the final iterative value, 


I” = calculated radiance for channel i from eq. (7) 
at iteration number n. . 


Then, following Smith [1970], in each channel the difference 
pnt a 
th (K) - BE(K)} (12) 


is independent of pressure within all atmospheric layers; 
hence the following iterative equation is obtained from 


eq. (11): 
=ntl =n n 


From a first guess temperature profile T(J), J = 1, 2, 
2.55 3,..., 100, Planck radiance values, B, (T(J) ] for all 
channels at each level can be computed from eq. (6). These 
values can then be used to determine layer mean Planck 
values, B, (K), by eqs. (8) and (9), and subsequently 


* radiances, I, > in accordance with eq. (7). The difference 
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between observed and calculated radiances, (ta, can 


be applied as a residual correction to obtain adjusted 
layer mean Planck values, BL) by use of eq. (13), and 
then improved radiances, ae can be calculated from 

eq. (7). 

By continuing the process of adjusting layer mean Planck 
values and calculating new radiance values until the differ- 
ence between observed and calculated radiance for each 
channel satisfies a convergence criterion at the final 
iteration step N, final layer mean Planck values, Bi(K) are 
obtained for all channels, i = 1,...,6 and all layers K = 
Legace < 

From the final layer mean Planck values a single mean 
temperature for each layer can be cinputed by noting that the 


fraction of calculated radiance for channel i and layer K, 


AI, (K), can be expressed 
N _ =N 
AT, (K) = By (K)4t, (K) ‘: (14) 


A similar expression applies to each of the six channels for 
a given layer, and a weighted layer mean Planck value can be 
formed by summing the at (K) over all i = 1,...,6. Therefore, 
for a layer K the weighted layer mean Planck value, By g(K) 
corresponding to a reference wave number, Vio for the layer 

is computed [Martin, 1973]: 


6 , 6 
=N _ —N 7 
Batak) = By (K) At, (K) . At, (K). | (15) 
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Substituting the weighted layer -mean Planck value into 


the Planck equation yields 


Byea(&) = o3,3/te RN -7 (16) 
where Ve is the reference wave number and T(K) is the mean 
temperature for layer K at the first guess step of iteration. 
The reference wave number for a layer K can be deter- 
mined by evaluating the weighted layer mean Planck value 
from eq. (16) using the initial layer mean Planck values 
Bi) (x), substituting the value of BO (x) thus obtained 
into a rearranged form of (16) 


(Co 3y/T(K)) 


542) (x) Le 


wid 1J - Cov 


2-0, (17) 
and solving for Ve using the Bailey iteration method of 
solution of transcendental differential equations [McCalla, 
1967]. Repetition for each layer gives reference wave 
numbers for all layers, Vy K = 1,2,...,17. 

The layer mean temperatures for each layer can then be 
calculated from the weighted layer mean Planck value corre- 


sponding to the final layer mean Planck value, By(K), by 


wtd 7] (18) 


where Vx is the reference wave number computed only at the 


initial step in accordance with (17). 
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Layer-mean temperatures may then be converted to layer- 


thickness values by integrating the hydrostatic equation 
between top and bottom pressure levels of the layer to 


obtain the hypsometric equation [Haltiner and Martin, 1957] 


AZ = Fa TK) in(-4) (19) 
& Po 
where 
AZ = thickness of a layer in meters, 
Ry.= 0.287 joules/(gm °K), 
g@ = 9.80 m/sec®, 


T(K) = mean temperature of the layer in °K, 


Pj» Po = pressures at top and bottomof layer in mb. 


B. APPLICATION 
1. Retrieval Input Parameters 

In addition to the satellite radiance data already 
discussed in Section III, first guess temperature profiles 
for each scan spot and atmospheric transmittance values for 
all six channels are required as input parameters for the 
retrieval program. 

a. First Guess Temperature Profiles 

A first guess temperature profile for each scan 

spot was derived from 56-level climatological profiles 
drawn from the U.S. Standard Atmosphere Supplement [1966]. 
The 15 N annual profile was assumed to be representative of 


an April profile at that latitude. Profiles for 30 N and 


30 





45 N for both January and July were interpolated with 


respect to time to give equivalent April profiles at each 
latitude. The three resulting "April" climatological 
profiles were then expanded to 100 level profiles by inter- 
polating temperature with respect to pressure to give 


2/7 or J-levels defined by eq. (4). 


temperatures at the p 

The 56-level climatological pressure levels and 
corresponding temperatures for the 15 N as well as the 
30 N and 45 N January and July climatological profiles are 
included in Appendix A, along with the interpolation. scheme 
used for expansion from 56 to 100 level profiles. 

Given the latitude of a scan spot, the corre- 
sponding first guess profile was obtained by interpolating 
with respect to latitude only between the 100 level "April" 
climatological profiles north and south of the scan spot, 
and by interpolating with respect to pressure between J = 2 
and J = 3 to obtain the J = 2.5 level temperature. The 
1000 mb climatological temperature was then replaced by the 
sea-surface temperature to "tie down" the profile at the 
lower boundary J = 100. 

b. Atmospheric Transmittance Values 
The absorbing gas structure for the six carbon 
dioxide channels was assumed to be represented by the 100 
J-level transmittances listed in Appendix B. These are 
transmittances calculated for a model atmosphere and the 


standard temperature profile which is also included in the 
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Appendix. However, as will be discussed in Section V, some 


of the transmittance profiles were later adjusted to improve 
retrieval results. 
2. Computational Procedure 
In practice, the reduction of VTPR "clear-column" 
radiances to atmospheric thicknesses was accomplished using 
two separate computer programs: one to retrieve mean 
temperatures for all scan spots considered, and a second to 
convert mean temperatures to thicknesses as well as to sort 
data by latitude band for analysis by the BIMED 02R regres- 
sion program. The two programs just described appear with 
sample outputs immediately following Appendix B. 
a. Mean Temperature Retrieval 

The procedure for retrieving layer mean tempera- 
tures for layers K = 1,2,...,17 from the VIPR "clear column" 
earbon dioxide radiances can be summarized as follows: 

[1] Derive a first guess temperature profile 
T(J), J = 1,2,2.5,3,...,100 in the manner described 
previously, and compute Planck radiance values, B, (T(J) ], 
for each J-level using eq. (6). 

[2] Calculate layer mean Planck values, 
BO) (x), K = 1,2,...,17 in accordance with eqs. (8) and (9). 

[3] Compute reference wave numbers for each 
layer using the initial layer~mean Planck values to form 
the weighted layer mean Planck value, BA cK) by use of 
eq. (15) and then solving (17) for Vx by the Bailey 


iterative method. 
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[4] Use the layer-mean Planck values to calcu-~ 


late the nth iterative radiances, Tye i =1,2,...,6 by eq. 
(7). | 

[5] Compare the observed radiances with the 
calculated radiances and apply the difference, [I, -Iy], as 
a residual correction to adjust the layer-mean Planck 
values in accordance with (13). . 

[6] Repeat steps [4] and [5] until convergence 
is achieved, convergence being defined as that condition in 
which 


N 


i < 0.0001 (20) 








1 


[7] From the final layer-mean Planck values, 
Bi (K), compute the corresponding layer-mean temperatures 
using eq. (18). 

b. Standard-layer Thickness Calculation 

Instead of calculating thicknesses for the 17 
layers for which mean temperatures were retrieved, it was 
decided to compute thicknesses between standard pressure 
levels that would permit a more complete and orderly exami- 
nation of layer combinations for specification of the 
troposphere and stratosphere as will be discussed in Section 
VII. The standard pressure levels chosen were the 
56-climatological profile levels plus the additional four 


levels at 650, 550, 450, and 80 mb. 
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Computation of thicknesses for the larger 


number of pressure intervals was accomplished by the 
procedure summarized as follows: 

[1] Consider the 17-layer mean temperatures 
obtained by retrieval eo located at the mid-levels of 
each layer, that is at levels J = 2.5, 7, 13, «s+, 97. 

[2] “Tie down" the resulting profile to the 
sea surface temperature at level J = 100 (1000 mb). 


[3] Determine the lapse rate between levels 


Cy 
" 


7 and J = 2.5 and continue the lapse rate to level 


01 (.01 mb) to determine the temperature at the top of 


Cy 
" 


the atmosphere. This completes a 19-level profile with 
temperatures at J-levels 1, 2.5, 7, 13, ...,; 97, 100. 

[4] Interpolate for temperature with respect to 
pressure, to derive temperatures at the 60 standard pressure 
levels from the 19-level profile on the J-scale. 

[5] Caleulate thicknesses of layers between 
successive standard pressure levels using a modification of 


the hypsometric equation (19) 


Re (T+ Ts) p ; 
ine. i CD (21) 
& 2 Po 


where 


ly T= temperatures (°K) at top and bottom 


Le “2 
pressure levels of the layer under consideration. 
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IV. EVALUATION OF RETRIEVAL ACCURACY 


If matching radiosonde soundings had been available, the 
retrieved mean temperatures foragiven "clear-column" scan 
spot could have been compared to mean temperatures calcu- 
lated from the corresponding radiosonde profile. However, 
matched radiosonde data were not available. 

Therefore to evaluate accuracy of the retrieved mean 
temperatures, profiles constructed from the retrieved 
temperatures Were first compared with the climatologically 
derived first-guess profiles to determine whether or not 
there were systematic differences between the two. Next, 
lapse rates between the bottom and top of layers of thickness 
AJ = 3 were examined and layers having super-adiabatic lapse 
rates were identified. 

Retrieved profiles were constructed by assigning mean 
temperatures to pressure levels at the middle of each layer 
by extending the lapse rate between levels J = 7 and J = 2.5 
to level J = 1 to compute a temperature for the top of the 
atmosphere, and by "tying down" the profiles by setting the 
temperature at J = 100 equal to the sea surface temperature. 
With temperatures then fixed for the 19 levels J = 1, 2.5, 
7, 135, «s«5 975 100, the temperatures at convenient inter- 
mediate levels were computed by interpolation to give the 
profiles T(J); 0 = 1, 2.5, 45 Ty aceg G4, 975 100, defined 


now at 34 J-levels for each scan spot. 
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Initial retrieval attempts resulted in profiles that 


were systematically colder than first guess profiles by a 
few degrees at all levels except J = 100 where the retrieval 
temperature was set equal to the sea surface temperature. 
The important bias upon retrieval was that super-adiabatic 
lapse rates of from 12 to 18 °K per km were observed between 
1000 and 900 mbs in virtually every retrieved temperature 
profile. Lapse rates above 900 mbs were approximately the 
Same as those of the profiles derived from climatological 
standards. | 

Similar systematic temperature errors were observed in 
retrieval methods using SIRS-A data [Fritz, Wark, et al., 
1972] and VIPR data [Jastrow and Halem, 1973]. In both 
references, the errors were compensated for by adjusting, 
or "tuning", the transmittance values channel by channel 
until the differences between retrieved profiles and the 
verifying radiosonde profiles at island stations used as 
check-profiles were brought within acceptable limits. 
However in this study, without radiosondes to use as check- 
profiles it was not possible to attempt "tuning" in a 
sophisticated manner. 

Nevertheless, to maintain quality control in the 
retrieved profiles it was necessary to eliminate the super- 
adiabatic lapse rates between 1000 and 900 mbs. This could 
have been done more elegantly, however reducing trans-~ 
mittance values for channel five by a factor of .95 and 


values for channel six by a factor of .90 proved to be 
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sufficient for the present purpose. Typical examples of 


retrieved profiles before and after the "tuning" of the 

’ transmittances as just described are depicted in Figs. 6 
and 7. Note that the super-adiabatic lapse rates between 
1000 and 900 mbs have been corrected by the “Guninge”" 


process which will be discussed further in Section V. 
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FIG. 6. Low-latitude temperature retrievals using both 
untuned and tuned transmittances. Note elimination of 
super~adiabatic lapse rate in lowest 100 mb by the tuned 
transmittances. : 
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FIG. 7. Mid-latitude temperature retrievals using both 
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untuned and tuned transmittances. Note elimination of super- 
adiabatic lapse rate in lowest 100 mb by the tuned transmittances. 





V. TRANSMITTANCE TUNING 


Although there are other possible causes of systematic 
error in radiance retrieval results, such as improper 
instrument calibration, the largest single source of such 
errors is presumably due to uncertainties in the trans- 
mittance functions [Drayson, 1971]. 

Transmittances are computed from theoretical models of 
absorption band structure, and are uncertain by at least a 
few percent in each channel [Jastrow and Halen, 1973). In 
the case of transmittance functions for the 15 um carbon 
dioxide channels, much of this uncertainty is due to lack 
of precise knowledge of the intensities, half-widths, and 
shapes of the absorption lines of all molecular absorbers 
contributing to tT, (p). Such knowledge should also be 
available for the 14 um band of ozone, and for the 20 um 
pure rotational band of water vapor [Drayson, 1971]. The 
ozone band absorbs weakly in the upper stratosphere, while 
a weakly absorbing edge of the water vapor band is effective 
in the lower troposphere. 

In addition, the transmittance functions are known to be 
weakly temperature dependent. For a given standard tempera- 
ture profile the total transmittance of the atmosphere for 
VIPR carbon dioxide channels can be considered to be the 
product of the individual transmittance of carbon dioxide, 


ozone, and water vapor [McMillin, Wark, et al., 1973]: 


(p) = (p) + (p) > (p) (22) 
TAP "CO, P "04 P "HO P 
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“Since the transmittance values used in this study were 


carbon dioxide transmittances for a mid-latitude standard 
atmosphere profile and not subjected to the correction of 
(22), it was assumed that the standard listings of t(p) 
were not completely descriptive of the absorption profile 
for the real atmosphere on the day the VTPR measurements 
were made. These factors are justification for periodic 
tuning in general [Jastrow and Halem, 1973], and in 
particular were considered to be the source of the systematic 
negative temperature errors in the initial retrieved pro- 
files already discussed. Based on this assumption, the 
transmittances were "tuned" to correct the super-adiabatic 
temperature lapse rates between 1000 and 900 mbs. 

Since the atmosphere transmits more strongly in certain 
pressure intervals than others for a given channel, the 
temperature of a given layer can be adjusted by tuning the 
transmittances for channels which receive greater radiance 
contributions from the layer. The selective transmittance 
of the atmosphere is illustrated in Fig. 8 for the "untuned" 
transmittances. It is clear from the figure that in the 
layer 1000 to 900 mb the largest atmospheric transmittance 
is in channels five and six. 

From the quadrature form of the radiative transfer 
equation (7), for a given channel radiance value, I,, it 
follows that selectively decreasing the layer weighting 
values, At, (K) requires that the value of the layer mean 


Planck values, By(K), be increased. An increase in the 


41 






(mb) 


PRESSURE 


d TAU/ dJ 


FIG. 8. Weighting functions for 
VTIPR carbon dioxide channels. 


layer-mean Planck value will give a higher value for the 
weighted layer mean Planck for a specific layer evaluated 
by eq. (15), and that in turn will give a higher layer-mean 
temperature in accordance with eq. (18). 

Increasing the retrieved mean temperature for any layer, 
K, could therefore be accomplished by decreasing the At; (K) 


for channels receiving radiance contributions from the layer. 
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Since channels five and six most significantly affected 


the layer 1000 to 900 mbs (as indicated in Fig. 8), increasing 
the retrieved temperature of this layer by decreasing At, (K) 
and At, (K) was feasible. 

Rather than decrease transmittance values for channels 
five and six only in the lowermost layers, a constant 
fractional decrease of transmittances was applied at all 
J-levels thereby decreasing At. (K) and At 6 (K) for all 
layers, K. Decreasing the transmittances only in the lower- 
most layers would have increased the temperature of the 
layer 1000 to 900 mbs, but it also would nate had the effect 


of changing the shape of the ord) weighting functions for 





channels five and six. Since a main objective of this 
thesis was to determine which atmospheric layer thicknesses 
were best specified by the VIPR "clear column" radiances, a 


tuning method which did not change the layers in which the 


at(J) 
ad 
absence of radiosonde check~-profiles to justify more complex 


curves "peaked" was considered desirable, in the 


tuning. The fractional decrease of transmittances at all 
levels satisfied this condition and had the effect of 
shifting the arts) curves of Fig. 8 for channels five and six 
slightly to the left while retaining their original shape. 

The selection of transmittance tuning factors, .95 for 
channel five and .90 for channel six, was strictly empirical. 
However subsequent to completion of the retrieval computations, 
it was learned through private communication [Dr. L. M. 


McMillin, 1973] that NESS had made somewhat similar tuning 


adjustments for these two channels for the same VIPR data 


analysis. 
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VI. RETRIEVED THICKNESS ANALYSIS > 


A. INDIVIDUAL LAYER SPECIFICATION 


Individual layer thicknesses were subjected to stepwise 
multiple regression analysis to determine which particular 


atmospheric layer thicknesses were best specified by using 





the “clear-column" radiance measurements as predictors. 


Prior to statistical analysis, retrieved thicknesses 
were separated into three samples based on the latitudes of 
their corresponding VTPR scan spots. 


number of scan spots for which thicknesses were retrieved 


are listed in Table 2. 


TABLE 2. 


Latitude bands and 


Sean spot samples for statistical analysis of 
retrieved thickness values. 


Latitude Band Latitude of Scan Spots 


La 


[Dixon, 1966] was used to analyze thickness data from the 


15 N 


JA 


Lat < 25 N 


25 N 


[A 


Lat < 35 N 


35 N Lat < 45 N 


JA 


Stepwise Regression Analysis 


The Stepwise Regression Analysis Program BIMED 02R © 


three latitude band samples. 


4y 


Sample Size 


82 


104 


LO9 





BIMED G2R computes, in a stepwise manner, a sequence 


of linear regression equations, with one variable added to 
the regression equation at each step. The variable added 

is the one that results in the greatest reduction in the 
previously unexplained sum of squares. This is also the 
variable which has the highest partial correlation with the 
dependent variable at the particular step in the analysis of 
variance. In addition, it is the variable which would have 


the highest F,-value when added at step k. 


k 


The FL -value at each step k is [Crow, et al., 1955] 


é 


where 


%(C.E.V.k) is the percent cumulative explained variance 
at step k. 


%(C.E.V.k-1) is the percent cumulative explained variance, 
step k-l. 


%Z(U.E.V.k) is the percent unexplained variance at step k. 


This study employed a statistical model expressed as 


sg = CC, + CN, t+ CaNe fF CaNs + CyNy + CON, + CeNe + CoNa 


0 a2 ca 
(25) 


where Az is a standard layer thickness, Ny through Ng are 
radiance predictors corresponding to VIPR "clear-column" 


radiance measurements of channels one through six, and 
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is that computed for the window channel (channel eight) 


Ny 
from the sea surface temperature by assuming a transmittance 
of unity [McMillin, Wark, et al., 1973]. Cy through Co are 
regression coefficients computed by BIMED O2R. 

For examination of the troposphere, the dependent 
variables, or predictands, were all possible layer thick- 
nesses between the pressure levels considered from 1000 to 
100 mb. For examination of the lower stratosphere, 
predictands were all possible layer thicknesses from 100 
to 20 mb. Each dependent variable was subjected to the 
BIMED 02R stepwise multiple regression analysis to determine 
how well the layer thickness was specified by the independent 
variables, the "clear~column" radiances. 


Related statistical parameters from the BIMED 02R 


output included: 


a. multiple R 

b. standard error of estimate, S.E. 
e. mean value, Az 

ad. standard deviation, a 

e. F.-value 

f. Re 


The most significant parameters for purposes of this study 


2 


were R® and S.E., which are related as follows: 


(S.E.)* = o@[(n-1)/(n-k-1) ](1-R?) (26) 


. where 
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n 9 : 
ao“ = z (Az, - Az)“/(n-1) (27) 


is the variance of a layer thickness for the sample, and 
n = sample size 
i = sample-element identifier 
k = number of predictors selected (k = 7 in this study). 
The fractional unexplained variance (LoR), can be 


approximated by 
(1-R?) £ GE)? (28) 


since (n-k-1)/(n-1) is close to unity for the sample 
sizes and number of predictors considered. 
The fractional explained variance can then be 


expressed 
eg = CS (29) 


which becomes the percent explained variance upon multi- 
plying R° by 100. 

The greater the percent explained variance for a 
particular layer thickness, the better the thickness is 
considered tobe specified by the radiance predictors. 

2, Individual Layer Results 
Although the specification of individual layers 


varied somewhat from latitude band to latitude band, in 
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general, the layer thicknesses best specified by the VTPR 


"clear-column" radiances were for relatively large pressure 
intervals, where also the lower pressure level was in the 
lower troposphere. Layer thicknesses most poorly specified 
were for pressure intervals in vicinity of the tropopause. 

A total of 99 different layer thicknesses were 
examined between 1000 and 20 mb for each latitude band. 
Rather than tabulate the statistical parameters for all 99 
layers, a representative group of layer thickness results 
has been extracted which illustrates the general trends in 
thickness specifications. 

Statistical parameters for a set of sequential 
layers between 1000 and 20 mb for all three latitude bands 
are listed in Table 3. Note that layers in vicinity of the 
tropopause, that is those between 150 and 90 mb, were not 
as well specified in terms of fractional explained variance 
as the layers above and below. This is indicative of 
"noise" in the retrieved mean temperatures near the tropo- 
pause associated with the reversal of temperature gradient 
across the interface. Both Tables 3 and 4 show that the 
explained variance of the thickness of layers in vicinity 
of the tropopause (e.g. 150-100, 100-90 mb) increases with 
increased latitude. Statistically this was a result of the 
larger standard deviations in the thicknesses of layers 150 
to 100 mb and 100 to 90 mb as one progresses into mid- 
“latitude bands. Give a larger 0°, the variable Az is more 
predictable in terms of layer-mean temperature and/or 


radiances. 
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TABLE 3. Statistical parameters of some representative 
sequential layers. 


150-100 
200-150 
250-200 
300-250 
400-300 
500-400 
600-500 
700-600 


1000-700 


Band 1 
S.E. R 
(gpm) 

1.670 -9TT9 

2.123 —-.9810 

6.175. — .8110 

7.966 .6627 

3.792 6052 

9.352 6076 

1.594 9317 

0.489 9944 

0.748 .9863 

1.572 9803 

1.255 .9821 

0.957 .9861 

2.532 -9T1T 

2.532 9797 


Band 2 
cao R° 
(gpm) 

1.487 9856 

2.981 9693 

5.731 8622 
6.222 8083 
2.914 7847 
8.907 7487 
3.686 8021 
0.844 9876 
0.514 9956 
1.436 9897 
1.469 .9853- 
1.415 9816 
1.403 9769 
2.929 .9813 


Lie) 


Band 3 
S.E. Re 
(gpm) 

1.661  — .9865 

2.940  .9758 

5.474 8982 

5.650 .8662 

2.623 8559 

8.558 .8622 

3.752 9283 

0.802  .9954 

0.532 9978 

1.504 .99h5 

1.531  .9921 

1.487  .9901 

1.467 9877 

2.960  .9894 





TABLE 4. Statistical parameters of some layers of special 


2 interest 
Band 1 Band 2 Band 3 
Pressure 2 2 2 
Layer S.E. R S.E. R S.E. R 
(mb) (gpm) (gpm) (gpm) 


300-150 0.384 = .9995 4.053 .9638. 4.072 9862 
400-150 1.230 .9982 2.826 9932 2.724 9974 
500-150 2.496 9957 1.778 .9985 . 1.594 9995 
600-150 3.452 9944 1.505  —.9993 1.462 .9997 
700-150 4.540 .9927 2.338  —.9987 2.446 .9994 


150-100 9.352  .6076 8.907  .7487 8.558 8622 
- ; 200-100 11.006 6565 12.603 =. 7383 12.306 .8813 
250-100 10.520 =. 7481 13.319 —. 7679 13.084  § .9105 
“100-90 3.792 6052 2.914 7847 2.623 8559 
100-70 11.704 6449 9.136 8008 8.274 8628 


1000-300 7.399 9808 8.668  .9833 8.939  .9909 
1000-150 7.090° .9989 5.011 9892 5.156 9982 
1000-100 2.390 .9988 4,351 =.9972 4.2u5 .9989 
150-20 27.430 8435 26.164 .8787 24.616 9144 
100-20 17.445 .9024 16.230 = .9239 15.0312 .9446 








. The sequential layers best specified were those 


between 600 and 200 mb. These layers are below the highly 
variable tropopause region, and furthermore are in a pressure 
interval where significant transmittance of radiance in 
channels three through six occurs (see Fig. 8). Presumably, 
having significant radiance contributions in a, given layer 
by several VTPR channels improves specification of the 
layer thickness by the radiance measurements. 

Statistical parameters for thicknesses of layers 
of special interest are listed in Table 4. These include 
the layers generally best specified, the evens generally 
most poorly specified, the 1000 to 300 mb layer thickness 
currently used in the FNWC vertical-analysis scheme [Holl ° 
et al., 1964], and the defined "tropospheres" and "strato- 
spheres" to be discussed in Sections VI (B) and (C). As 
in the case of sequential layers already discussed, the 
most poorly specified of all pressure intervals in Table 4 
were shallow and had boundaries near the tropopause. The 
best specified layers were those in the mid-troposphere where 
Significant radiance sonteipurtons from at least four of 


the co, channels occur. 


B. TROPOSPHERIC SPECIFICATION BY SUB-LAYERS 

Results of the statistical regression analysis applied 
to individual layers were next used to examine various com- 
binations of the sequential sub-layers in the troposphere to 
determine which of the combinations could best specify the 


full tropospheric thickness. 
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1. Method of Analysis 


’ To compare the results of different tropospheric 
thickness combinations, a normalized root mean square error 
index was defined. The error index, E.I., for a two-layer 


combination, for example, was 


E.I. = : (30) 
2trop 





where S.5.. and S.Bs, are the standard errors of estimate 


1 2 


for the individual layers, and Bes is the mean thickness 
of the defined troposphere. Similarly, the E.I. for a 


three-layer combination was 


2 
E.I. = i eae eee (31) 
2trop 


Weighting of standard errors when computing the E.I. for 
a particular layer combination was not required since the 
standard error of thickness reflects variance in the retrieved 
mean temperatures, and weighting of mean temperatures by 
In(p,/P5) factors was implicit in conversion of mean 
temperatures to thicknesses by the hypsometric equation. 

This error index expresses combined layer errors 
in a root mean square (RMS) sense, and the smaller the E.I. 
value the better the combination of layers specifies the 


total thickness. 
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Specification of the tropospheric thickness was 


considered by defining the troposphere in one series of 
tests to be the layer 1000 to 100 mb, and in a second series 
of tests to be the layer 1000 to 150 mb. In each case, best 
specification of the "troposphere" was sought by comparing 
error-index values for various sub-layer combinations within 
the troposphere. 

To form combinations of layers, the concept of a 
"sliding layer" was adopted. The "sliding layer" was first 
tied down tothe 1000 mb level-and the RMS error, as expressed 
by the error index, was computes for the resulting two-layer 
combination. The sliding layer was then stepped to a higher 
level and the RMS error for the resulting three-layer combina- 
tion calculated. The procedure of stepping the sliding 
layer to higher levels to form new combinations was continued 
until a final two-layer combination was formed with the top 
of the sliding layer being coincident with the defined tropo- 
pause level (150 or 100 mb). The pressure interval of the 
sliding layer was then increased and the procedure repeated. 
"Sliding layers" of 300, 400, 500, 600, 700, and 800 mb 
were employed. Fig. 9 illustrates the combinations of sub- 
layers formed by the 400 mb sliding layer in the 1000 to 
100 mb troposphere analysis. 

Error index values from all combinations of layers 
formed by the above procedure were compared to determine 
. which particular combination gave the minimum RMS error for 


specification of the full tropospheric thickness. In 
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400 






PRESSURE (MB) 


1 2 3 4 5 6 


CASES 


FIG. 9. All combinations, or cases, of tropospheric 
sub-layers formed by the 400 mb "sliding layer" with 
the tropopause defined at 100 mb. In each case, the 
combined sub-layers span the pressure interval 1000- 
100 mb and the location of the sliding layer is 
indicated by the solid line. 
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addition, the combination layer results were compared 


against the E.I. values for the "tropospheres" considered 
as single layers, 1000 to 150 mb and 1000 to 100 mb. 

2. Troposphere Results 

| In all three latitude bands, defining the tropo- 
sphere to be the pressure interval 1000 to 150 mb, rather 
than 1000 to 100 mb, resulted in generally lower E.I. values 
for the layer combinations considered. This result can be 
attributed to the fact that the single layer 150 to 100 mb 
was the most poorly specified of all troposphere layers due 
to variability induced by the tropopause in this layer. 
Therefore defining the model tropopause at 150 mb gave 
smaller combined error index values by eliminating errors 
of the 150 to 100 mb layer. 

In latitude bands two and three, the maximum RMS 
error for both two and three layer combinations occurred 
when the lower limit of the top layer of a combination was 
at 250 mb. In latitude band one this maximum Secured when 
the lower limit of the top layer was at 200 mb. This result 
may also be explained in terms of influence of the tropopause, 
since the tropopause is generally higher in lower latitudes 
[Haltiner and Martin, 1957], and therefore its detrimental 
effect on the standard error of estimation is greater for 
the layer 200 to 100 mb than for the layer 250 to 100 mb 
as in the higher latitude bands. Standard error values for 


these layers listed in Table 4 illustrate this relationship. 


20 


The general trend in the troposphere was for the 
combination error index values to be lowest for certain 
two-layer combinations topped at 150 mb with a dividing 
level near 650 mb (Table 5). Table 5 lists in sequence the 
layer-combinations found here to be most effective in giving 
minimum error index values. Examination of Figs. 10, 1l, 

12 show that other layer choices give rise to larger RMS 
errors of specification in the troposphere. No attempt has 
been made in this study to ascribe statistical significance 
to any one layer choice because of the limited period of 

the test. However, the general nature of the "best" results 
in each band seem to resemble the tropospheric stratification 
already described. 

When the RMS errors for layer combinations were 
compared against those for the "tropospheres" taken as 
single layers, it was found that the layer 1000 to 150 mb 
could be better specified by certain combinations of sub- 
layers than when considered as a single layer. ‘The layer 
1000 to 100 mb, however, had the minimum of all RMS errors 
found when considered as a single layer. This evidently 
was due to the fact that the relatively poor specification 
in the 150 to 100 mb layer became insignificant in compari- 
son to the thickness-specification of the layer 1000 to 
100 mb. Error index values for the two defined troposphere 
layers considered as single layers are listed by latitude 
band in Table 6, and these may be compared with the best 


combination results in Table 5. Note that the best 
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TABLE 5. Tropospheric sub-layer combinations having the 
smallest error index values. 


Pressure Layers (mb) Band 1 E.I. 


1000-700, 700-400, 400-150 .000315 
1000-700, 700-100 . 000341 
1000-550, 550-150 000370 
.Band 2 
1000-700, 700-150 | .000275 
1000-650, 650-150 000308 
1000-600, 600-150 .000338 
| Band 3 
1000-700, 700-150 000284 
1000-650, 650-150 . .000316 
1000-600, 600-150 000345 


oy | 





specification of the troposphere in each latitude band 


was obtained by considering the single alyer 1000 to 100 mb. 


TABLE 6. Error index values for the defined "tropospheres" 
considered as single layers. 


"Troposphere" Band 1 Band 2 Band 3 
(mb) E.I. ELI. E.I. 

1000 to 100 .000148 000271 .000265 

1000 to 150 -000515 .000364 .000382 


Error index values for all tropospheric layer 
combinations considered, as well as values for the defined 
tropospheres taken as Single layers, are graphically 
depicted in Figs. 10, 11, and 12. For each latitude band, 
graphs (a), (b), (c) depict E.I. values for layer combina- 
tions formed by "sliding layers" of 300, 400, and 500 mb. 
The (d) graph for each band depicts E.I. values for all 
other combinations considered, that is combinations 
resulting from sliding layers of 600, 700, and 800 mb, as 
well as error index values for the single layer troposphere 
cases. 

Note that except for the two-layer combination 
1000 to 700 and 700 to 100 mb in latitude band one, all 
other combinations formed by sliding layers had smaller 
error index values when the tropopause was defined at 


150 mb. Also note the shift in maximum E.I. value from 
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FIG. 10. Cases of tropospheric sub-layer combinations considered and corresponding error 
index values (15 N to 25 N). 

® E.I. values with tropopause at 150 mb. OE.I. values with tropopause at 100 mb. 
(a) Sub-layer combinations formed by sliding layer of 300 mb, (b) 400 mb, 

(c) 500 mb, (ad) all others. 
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1000 
FIG. 11. Cases of tropospheric sub-layer combinations considered and corresponding error 
index values (25 N to 35 N). 


# E.I. values with tropopause at 150 mb. 4 E.I. values with tropopause at 100 mb. 
(a) Sub-layer combinations formed by sliding layer of 300 mb, (b) 400 mb, 
(ec) 500 mb, (d) all others. 
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FIG. 12. Cases of tropospheric sub-la 
index values (35 N to 45 N). 
# E.I. values with tropopause at 150 mb. 4 E.I. values with tropopause at 100 mb. 


(a) Sub-layer combinations formed by sliding layer of 300 mb, (b) 400 mb, 
(c) 500 mb, (d) all others 
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combinations with the topmost layer having a lower limit 


at 200 mb in band one to combinations with the topmost 
layer having a lower limit at 250 mb in bands two and 


three. 


C. STRATOSPHERIC SPECIFICATION BY SUB-LAYERS 
Best specification of the stratosphere to 20 mb was 
sought in a manner similar to that of the troposphere 
already discussed. 
1. Method of Analysis 
The error index was as defined in specification of 
the troposphere, except that Az was used instead of 


strat 
in eqs. (30), (31). The stratosphere was considered 


AZ nop 
in one series of tests to be the layer 150 to 20 mb, and in 


/ 8 Second series of tests to be the layer 100 to 20 mb. 
"Sliding layers"of 30, 40, 50, 60, 70, and 80 mb were used 
to form layer combinations. Figure 13 illustrates the 
layer combinations formed by the 30 mb sliding layer 
between 150 and 20 mb. 
2. stratosphere Results 

Although the results do not afford as sharp a 
specification as for the troposphere, a general result was 
that the 100 to 20 mb "stratosphere" had lower RMS errors 
for layer combinations than the 150 to 20 mb case. This 
again can be attributed to the influence of the tropopause 
in the layer 150 to 100 mb. By defining the tropopause to 
be at 100 mb the poorer specification of the layer 150 to 


100 mb was eliminated. 
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FIG. 13. All combinations, or cases, of 
stratospheric sub-layers formed by the 30 mb 
"sliding layer" with the tropopause defined at 
150 mb. In each case, the combined sub-layers 
span the pressure interval 150-20 mb and the 
location of the sliding layer is indicated by 
the solid line. 
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The three best layer combinations in each latitude 


band are listed sequentially in Table 7 with their corre- 
sponding error index values. There was not a unique set 
of combinations that proved to be most effective for all 
latitudes, but a general. trend for the error index of best 
combinations to decrease with increasing latitude is evident. 
This apparently was a result of the tropopause being at 
lower levels for higher latitudes, which tended to reduce 
the influence of the tropopause on specification of the 
defined stratosphere above 100 mb at higher latitudes. AS 
in the tropospheric analysis, no attempt has been made to 
attach statistical significance to the results due to the 
limited period of the test. 

Error index values for the two "stratospheres" 
considered as single layers were not as good as for 
combinations of layers. A decrease in E.I. values when 
the 150 to 100 mb layer was icornonetea into the total 
thickness was not observed as in the troposphere analysis. 
This difference was evidently due to the fact that the 
thickness-error induced by the 150 to 100 mb layer was 
relatively smaller in comparison to the specification of 
the full tropospheric thickness, but had greater influence 
on thickness-error of the tiestcesnees., Error index values 
for the two defined stratosphere layers are listed by 
latitude band in Table 8. 

Error index values for all stratospheric layer 


combinations considered, as well as values for the defined 
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TABLE 7. 


Pressure Layers (mb) 
100-90, 90-60, 60-20 
100-80, 80-20 

150-100, 100-70, 70-20 


100-90, 90-60, 60-20 
100-80, 80-50, 50-20 


100-70, 70-20 


100-70, 70-30, 30-20 
100-80, 80-40, 40-20 
100-90, 90-60, 60-20 
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Stratospheric sub-layer combinations having the 
smallest error index values. 


Band 1 E.I. 
.00123 
.00130 
.00131 
Band 2 
.00105 
.00112 
-00116 
Band 3 
.00094 
.00095 
.00096 





TABLE 8. Error index values for the defined "stratospheres" 
considered as single layers. 


"Stratosphere" Band 1 Band 2 Band 3 
(mb) E.I. E.I. E.L. 
100 to 20 .00172 .00162 .00148 
150 to 20 . .00220 .00214 .00198 


stratospheres taken as single layers, are graphically 
depicted in Figs. 14, 15, and 16. For each latitude band, 
graphs (a), (b), (c), (d) depict E.I. values for layer 
combinations formed by 30, 40, 50, and 60 mb "sliding 
layers"; (e) graph depicts error index values for the 
single-layer stratosphere cases in addition to those for 
the 70 and 80 mb sliding layer combinations. 

The graphical display of stratospheric error index 
values illustrates the fact that defining the tropopause to 
be at 100 mb gave better E.I. values than when the tropo- 
pause was defined at 150 mb and the layer 150 to 100 mb was 
included in the stratospheric specification. Also note 
that the E.I. values are an order of magnitude greater than 
the tropospheric error index values previously examined. 

Apparently, the significantly greater E.I. values 
in the stratosphere are due to standard errors of estimate, 
S.E., of layers in the stratosphere reflecting uncertainties 
. of explained variance for the tropospheric column below, as 


well as for the layers of the stratosphere. In the 
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troposphere, on the other hand, the layer combinations are 


"tied down" to the surface and tropospheric standard errors 
reflect only unexplained variance in the troposphere. 
Essentially, the relatively greater S.E. values of the 
stratosphere result in correspondingly greater E.I. values, 
although Az of the denominator is somewhat smaller than for 


the tropospheric’ case. 
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FIG. 14. Cases of stratospheric sub-layer combinations considered and corresponding error 
index values (15 N to 25 N). 

» E.I. values with tropopause at 150 mb. 4 E.I. values with tropopause at 100 mb. 
(a) Sub-layer combinations formed by sliding layer of 30 mb, 

(bo) 40 mb, (c) 50 mb, (d) 60 mb, (e) all others. 
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FIG. 15. Cases of stratospheric sub-layer combinations considered and corresponding error 
index values (25 N to 35 N). 

=» E.I. values with tropopause at 150 mb. 4 E.I. values with tropopause at 100 mb. 
(a) Sub-layer combinations formed by sliding layer of 30 mb, 

(b) 40 mb, (ec) 50 mb, (a) 60 mb, (e) all others. 
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FIG. 16. Cases of stratospheric sub-la 
index values (35 N to 45 N). 

» E.1I. values with tropopause at 150 
(a) Sub-layer combinations formed by sl 
(b) 40 mb, (c) 50 mb, (da) 60 mb, 





yer combinations considered and corresponding error 


mb. 4 E.I. values with tropopause at 100 mb. 
iding layer of 30 mb, 
(e) all others. 





VII. CONCLUSIONS 


This thesis has demonstrated that atmospheric thick- 
nesses can be retrieved from VIPR "clear-column" carbon 
dioxide radiances by an improved iterative technique for 
direct solution of the radiative transfer equation, and that 
systematic errors in the retrieval results can be corrected 
by selective adjustment of the transmittances. In addition, 
the retrieved thicknesses of certain standard layers appear 
to be better specified than others when the "clear column" 
radiances are used as regression predictors. Furthermore, 
it has been found in this case that specification of the 
full troposphere and stratosphere layers can in some 
instances be improved by decomposition of the layer into 
sequential pressure intervals. 

It should be noted that the study conducted here bears 
more heavily on the retrieval technique, and utilizes 
simulated atmospheric-layer thicknesses which have had 
"clear-column" radiances built into them. In order to make 
any final inferences as to choice of an optimum layer 
thickness, the procedure suggested here would have to be 
conducted over a period of time and retrieved thicknesses 
verified against thicknesses computed from corresponding 


radiosonde soundings. 
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The following procedure was used to expand standard 


56 level climatological temperature profiles to 100 level 
profiles on the J scale: 

1. Where pressures of J-levels were equal to pressures 
of standard levels, the temperature at the J-level was set 
equal to the temperature at the standard level. 

2. Where pressures of J-levels were not equal to 
pressures of standard levels, temperatures at the J-levels 


were computed by interpolation as follows: 


sr Trp ep(a) = Par (MeL) 
STD STD : 
where 
T(J) = desired temperature at the J-level 


P(J) = pressure at the J-level 


(M-1)= known temperature at nearest standard pressure 
level above the J-level 


Topp (M) = known temperature at nearest standard pressure 
level below the J-level 
Pgpp(M-1) 
= pressures at nearest standard levels above and 
Pomp (™) below the J-level. 
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